Chapter 1 

Modelling Starburst in Hn galaxies: From 
chemical to spectro-photometric evolutionary 
self-consistent models. 



Martm-Manjon, M.L., Molla, M., Diaz, A.I. & Terlevich, R. 



Abstract We have computed a series of realistic and self-consistent models that 
reproduce the properties of Hll galaxies. The emitted spectrum of Hll galaxies is 
reproduced by means of the photoionization code CLOUDY, using as ionizing spec- 
trum the spectral energy distribution of the modelled Hll galaxy, calculated using 
new and updated stellar population synthesis model (PopStar) This, in turn, is cal- 
culated according to a star formation history and a metallicity evolution given by a 
chemical evolution code. Our technique reproduces observed abundances, diagnos- 
tic diagrams, colours and equivalent width-colour relations for local Hll galaxies. 



1.1 The Model. 

Hll galaxies are characterized by strong and narrow emission lines and by a low 
metal content, but this does not necessary mean that these galaxies be young sys- 
tems. The current burst of star formation (SF) dominates the spectral energy dis- 
tribution (SED) even if previous stellar populations are present, making difficult to 
know the star formation history (SFH) of the galaxy. We have made a new grid 
of star-bursting models, based on |6), using simultaneously the whole information 
available for the galaxy sample: the ionized gas, which defines the present time state 
of the galaxy, and spectrophotometric parameters, related to its SFH. These mod- 
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els imply the combination of three different tools: a chemical evolution code, an 
evolutionary population synthesis code and a photoionization code, all computed in 
a self-consistent way, that is, using the same assumptions regarding stellar evolu- 
tion, model stellar atmospheres and nucleosynthesis, and a realistic age-metallicity 
relation. 

The model consist in a set of successive burst of SF in a region with a total mass 
of gas of 100 10 6 Mq. Each model is characterized by the initial efficiency, e, the 
percentage of available gas consumed to form stars, and the attenuation: the initial 
efficiency is attenuated by a constant factor, k, in the successive burst following the 
expression y/n = yb(&) 1 > where n is the number of the current burst. 

We present two models with k=0.2 and two different initial efficiencies: Model 
1, with £=50% and Model 2 with £=12%. 

With the chemical evolution code, based on J3J, we obtain for every 0.7 Myr 
step the abundances of 15 elements: H, He, C, O, N, Ne, Na, Al, Mg, Si, S, Ca, Ar, 
Ni, Fe, the star formation rate (SFR) and the corresponding age-metallicity relation, 
Z(t). 

A SED from the Single Stellar Populations (SSPs) library PopStar Q, has been 
assigned to each time step taking into account the SFH, £, f (t) and Z(t) obtained 
from the chemical evolution model. When more than one burst take place, the SED 
is the sum of the SEDs of every stellar generation convolved with the SFH. The final 
result is the total luminosity at each wavelength for each time step, (t): the SED 
corresponding to the whole stellar population, including the ionizing continuum of 
the last formed stellar generation. 

The gas is ionized by the massive stars belonging to the current burst of SF. 
We use CLOUDY Q in order to obtain emission lines. Each burst photoionization 
model is characterized by a radius R, calculated according to the mechanical en- 
ergy output of the massive star winds and SNel explosions, a constant gas density 
n#=100 cm~ 3 , an ionization parameter U = Q(H)/47rcn#R 2 , where Q(H) is the num- 
ber of ionizing photons, and the chemical abundances obtained from the chemical 
evolution code. 



1.2 Results. 

In fig ll.ll (left) we represent the SFR of Model 1 and Model 2, and the evolution 
of the oxygen abundance along 13.2 Gyr (right). The first burst is strong, while 
subsequent are less intense due to the strong attenuation and the decrease of the 
available gas to form stars. We see in fig. 11.11 (right) that the efficiencies of the 
Model 1 and Model 2 are the upper and lower limits, respectively, for the oxygen 
abundance in Hll galaxies. 

The properties of the ionized gas, showed as excitation diagnostic diagrams, can 
be seen in fig J 1 .21 Model 1 reproduces the more excited and more metallic galaxies, 
with high [OIII]/H|3 due to its high efficiency of SF, while Model 2 reproduces less 
metallic galaxies, with high [OIII]/H|3 and low [OIIJ/Hp ratios. 
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Fig. 1.1 SFR along 13.2 Gyr of Model 1 and Model 2 (left). Evolution of the oxygen abundance 
for both models (right). Black dashed lines corresponds to the observed limits for Hll galaxies (4l. 
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Fig. 1.2 Diagnostic diagrams for low efficiency model (left) and high efficiency model (right) 



The broad band colours have been computed including the contribution by the 
stronger emission lines. These are the colours which are readily observable though 
integrated photometry (fig ll.3t . 
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Fig. 1.3 Colour-colour diagram for Model 1 (right) and Model 2 (left). 
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Most Hll galaxies have EW(H^) values lower than 150 A, which indicates the 
presence of an old non-ionizing population. According to the models, these values 
can not be produced by just one generation of stars (a SSP), and it is necessary to 
invoke an underlying population, of at least 1.3 Gyr, to contribute to the continuum. 
EW(Hj3) vs. (U-V) along the evolution of the galaxy for the stellar bursts is shown 
in fig J 1 .41 The trend of the data is reproduced by our models and not by SSP's since 
in that case the colours are bluer than observed. In order to decrease EW(H^) and 
to obtain redder colours, a more metal-rich SSP might be selected (dotted line at the 
left panel), but such high abundance does not reproduce the observations (fig. ll.ll ). 

In summary, our models are able to produce redder colours with the contribu- 
tion of previous stellar generations and reproduce better the trend of Hll galaxies 
observations, as well as abundances, colours and emission lines. 
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Fig. 1.4 (left) EW(H^) vs (U-V), compared with two SSPs from Starburst99 0; a very low metal- 
licity SSP (black solid line) and a high metallicity SSP (dotted line). The contribution of the under- 
lying continuum make the colour to be redder and reproduce observational trend of Hll galaxies, 
as we see on the right: EW(H^) vs a pseudo-colour equivalent to U-V. 
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